Collective rearrangement at the onset of flow of a polycrystalline hexagonal columnar 

phase 
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Creep experiments on polycrystalline surfactant hexagonal columnar phases show a power law 
regime, followed by a drastic fluidization before reaching a final stationary flow. The scaling of the 
fiuidization time with the shear modulus of the sample and stress applied suggests that the onset 
of flow involves a bulk reorganization of the material. This is confirmed by X-ray scattering under 
stress coupled to in situ rheology experiments, which show a collective reorientation of all crystallites 
at the onset of flow. The analogy with the fracture of heterogeneous materials is discussed. 

PACS numbers: 82.70.-y, 83.50.-v, 62.20.Hg, 62.20.Fe 
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Yield stress fluids are materials that respond elastically 
below a yield stress, a y , and flow as liquids for stresses 
above a y Typical materials that fall in this category 
include granular media, dry foams or dense colloidal sus- 
pensions and emulsions. Very recently, Becu et al. [2[ 
have shown that the behavior of concentrated emulsions 
at the onset of flow depends crucially on the microscopic 
interactions between the droplets, thus underlining the 
non universal nature of the solid to fluid transition of 
yield stress fluids. Several fundamental issues are still 
unclear regarding this transition. These include in par- 
ticular the possible prediction of if and when a material 
will flow, from its behavior prior to flow. In addition, 
whether the solid to fluid transition is accompanied by 
structural modifications of the yield stress fluids remains 
largely unknown. 

We perform creep experiments in order to investigate 
the onset of flow of soft polycrystalline hexagonal colum- 
nar phases. Upon application of a small and constant 
stress, <7, the material deforms plastically in a similar 
manner as a solid By contrast, under high stress, 
it behaves as a fluid and flows with a constant shear 
rate that depends on a 0, Quite remarkably, for 
an intermediate value of er, an intriguing regime is ob- 
served during which the material, after an incubation 
time that depends on both a and the sample shear mod- 
ulus, evolves suddenly from a power law creep to a flow 
regime, through a dramatic decrease of the viscosity. 
Thanks to a combination of Synchrotron X-ray scattering 
experiments under stress and in situ rheology measure- 
ments, we are able to correlate the solid to fluid transition 
with structural modifications of the material, thus shed- 
ding light on the mechanisms involved at the onset of 
flow of some complex fluids. 

The experimental system is a lyotropic hexagonal 
phase consisting of infinitely long oil tubes which are sta- 
bilized by a surfactant monolayer and immersed in wa- 
ter. The tubes of uniform radius arrange on a triangular 
array and form, therefore, a two-dimensional columnar 
crystal. The samples' composition and elasticity have 



been described elsewhere 0, Q ■ Four samples with vari- 
able shear moduli, Go, ranging from 250 to 3200 Pa, have 
been investigated. Rheology experiments are performed 
on a stress-controlled Paar Physica UDS 200 rheometer 
in a Couette geometry with a gap of 1 mm. Once loaded 
in the cell, the sample is systematically cooled below the 
hexagonal to isotropic fluid phase transition temperature 
(~ 10° C), and then warmed up to room temperature. A 
polycrystal, with an isotropic distribution for the crys- 
tallites orientation, is thus rapidly recovered. At time 
t = 0, we then apply a constant shear stress, a, to the 
material and record the time evolution of the strain, 7. 
An instantaneous jump of 7 is measured, revealing the 
elastic nature of the material at high frequency, which is 
followed by a monotonic increase of 7. The evolution of 
the instantaneous shear rate, 7, calculated as the time 
derivative of 7, is shown in fig. [T] for several stresses, and 
displays three distinct regimes. Initially, 7 decreases as 
a power law with time, indicating a continuous increase 
of the instantaneous viscosity, n, defined as cr/7. This 
power law regime is followed by a very abrupt increase of 
7 of more than 2 orders of magnitude, until the shear rate 
reaches a constant and high value, which corresponds to 
a stationary flow with a lower viscosity. For t larger than 
0.1 sec (the time resolution of our experiments), and over 
at least three orders of magnitude, we find in the power 
law regime 7 ~ t p with p — —0.34 ±0.05, whatever a and 
Go- This regime recalls the primary creep regime, the 
so-called Andrade creep Q measured in many crystalline 
and non-crystalline materials. However, the catastrophic 
acceleration of the shear rate that follows the primary 
creep, which illustrates a drastic decrease of the viscos- 
ity, is neither observed for metallic materials, nor for 
yield stress fluids such as foams |8| or colloidal glasses 
[9(, where the transition from creep to flow is smooth. 
The features shown in fig. [1] appears rather unique for 
a complex fluid, but in fact display intriguing analogies 
with the fracture under constant load of many composite 
heterogeneous materials 13, U|. We define the fluidiza- 



tion time, Tf, as the time at which 7 is minimum (i.e. 
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FIG. 1: Time-evolution of the shear-rate in creep experiments 
for a sample with an elastic shear modulus Go = 2240 Pa. 
At time a constant stress a = 12, 13, 15, 19, 21, 22, 30, 35 Pa 
(from bottom to top) is applied to the sample. 



the instantaneous viscosity is maximum). 

We report the stress dependence of the fiuidization 
time in fig. [2^, for several samples with different elas- 
ticities. We find that, for all samples, r/ decreases as a 
increases, and that curves for softer samples are shifted 
to lower stresses. Interestingly, all data points can be put 
onto a single master curve if the stress is renormalized by 
a power law of the shear modulus, G™. We find the best 
collapse of the data for m = 0.64. For this value, the 
master curve can be reasonably well fitted with a power 
law T f ~ (cr/Gg- 64 )(- 31±0 - 2) (fig. [S)). This scaling im- 
plies that there is no divergence of Tf with stress. It 
indicates that the material does not show any detectable 
yield stress in our experimental time window, in contrast 
with the findings for micellar polycrystals for instance 
On the other hand, we have shown previously that 
the shear modulus of the polycrystal scales as R~ 3 , where 
R is the radius of the surfactant tubes Q. Hence, the nor- 
malized stress <t/Gq scales roughly as a x R 2 . We ex- 
pect the number density of tubes in the material to scale 
as 1/i? 2 . Thus, the master curve indicates that Tf de- 
pends only on the stress applied per tube. Interestingly, 
despite variation of Tf of two orders of magnitude, the 
strain 7/ at time t = Tf is approximately constant. We 
find: 7/ ~ (69 ± 23)%, irrespective of the stress applied 
and of the sample elasticity (fig. Eb). This shows that the 
onset of flow results from a strain controlled mechanism. 
Note that a critical strain similar to ours (i.e. of the or- 
der of 1) is typical for shear-induced order in hard sphere 
suspensions [13| or to completely fluidize a foam sample. 
Both a critical strain of 1 and the remarkable dependence 
of the stress with the number of tube suggests that the 
fiuidization time reflects a bulk mechanism that involves 
the whole sample and not an interfacial property. This 
is supported by the results described below. 

In order to gain a better understanding on the mecha- 
nisms at work in the creep regime and at the onset of flow, 




FIG. 2: Fiuidization time, Tf, as a function of (a) the stress 
applied a, and (b) a/Go' 64 (symbols), for samples with dif- 
ferent shear moduli, Go = 250 Pa (squares), 760 Pa (crosses), 
2240 Pa (triangles) and 3200 Pa (stars) ; the line in (b) is a 
power law fit yielding an exponent —3.1 ± 0.2. (c) Strain at 
t — Tf as a function of a/Go' 64 (same symbols as in (a,b)). 



we have coupled time-resolved Synchrotron small-angle X 
ray scattering (SAXS) under stress to in situ rheological 
measurements. Experiments have been performed on the 
ID-2 beamline at the ESRF, Grenoble, France. A stress- 
controlled Haake RS300 rheometer equipped with a Cou- 
ette cell is used. The scattering profiles are recorded for 
two geometries, the incident beam being either radial or 
tangential with respect to the cell. Only the softer sam- 
ple (Go = 250 Pa) has been investigated. A procedure 
analogous to the one used for the rheology experiments 
has been performed in order to ensure a reproducible 
polycrystalline initial state corresponding to an isotropic 
distribution for the orientation of the crystallites, as con- 
firmed by isotropic rings found in the SAXS patterns for 
both radial and tangential geometries. We can evaluate 
that the size of the crystallites is of the order of 1 fim 
[3]. Taking into account the size of the incident beam 
(100 /im) and that of the gap (1 mm), we conclude that 
the SAXS spectra reflects the orientation distribution of 
a large assembly of crystallites. 

A typical series of radial scattering patterns are shown 
in fig. [3^ together with the strain, 7, and the instanta- 
neous viscosity, 77. The time evolution of 77 clearly shows 
the three successive regimes (power law creep, fiuidiza- 
tion and stationary flow) described previously. In fig. [3k., 
the time is normalized by t c , the time at which the mate- 
rial starts to flow with a stationary shear rate. The SAXS 
patterns exhibit isotropic rings at the beginning of the ex- 
periment, and then converts into a markedly anisotropic 
signal in the flow regime, with higher intensity along the 
vorticity direction (6 = 90 and 270 deg) . To quantify the 
time evolution of the patterns and correlate it with that 
of the rheological properties, we compute the angular de- 
pendence of the scattered intensity integrated around the 
first Bragg peak (azimuthal scan), and identify the angu- 
lar position of the local maxima. Figure [3)3 displays the 
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FIG. 3: Combined rheological and structural measurements 
in a creep experiment. The sample elasticity is Go = 250 Pa 
and the stress is a — 2 Pa. The time is normalized by the time 
t c at which the sample flows with a constant shear rate. The 
dashed vertical line marks the beginning of a stationary flow 
regime, (a) In situ rheological measurements, strain (thick 
line) and instantaneous viscosity (thin line), and SAXS pat- 
terns taken from left to right at t/t c ~ 0.25,0.75,0.98 and 
1.14. The grey scale is linear and identical for all images, (b) 
Time evolution of the angular position of the local intensity 
maxima in an azimuthal scan. The bubble size scales with 
intensity, (c) Time evolution of the contrast C (line) and the 
scattered intensity at 9 = (180 ± 2) deg (solid circles). 



evolution of the position and intensity of these maxima, 
for 9 in the range 160 — 340 deg (9 = and 180 deg. corre- 
sponds to the velocity direction). The experimental data 
clearly show that, as soon as the systems flows (t/t c > 1), 
the surfactant tubes are in majority aligned along the 
flow (maxima around 270 deg) , although the angular dis- 
tribution for the tubes orientation is wide (A9 ~ 35 deg). 
We note that the angular width, A9, and the shear rate 
measured in the flow regime (7 = 0.14 sec -1 for a = 2 
Pa) are in quantitative agreement with our previous find- 
ings 0, [fl concerning both the continuous decrease of A9 
as 7 increases and the stress dependence of 7, suggest- 
ing that the flow regime investigated here belongs to the 
shear-thinning regime previously determined for 7 in the 
range (1 — 40) sec -1 . We calculate the contrast, C, de- 



fined as (I m ax ~ Irnin) / (I-max + I-min) i where I m ax and 

Imin are respectively the maximal and minimal intensity 
along the azimuthal scan profile ; C varies between 0, for 
a perfectly isotropic distribution of the crystallite orien- 
tation, and 1 for a totally anisotropic distribution. As 
shown in fig. [3J:, C is small (~ 0.25) at the beginning of 
the experiment, and after a slight decrease [lq . increases 
continuously, until it reaches a constant and high value, 
C ~ 0.92, at t = t c . 

Because each experiment probes a given ensemble of 
crystallites, the exact patterns of the peak position along 
a scan change from one experiment to another one ; in 
particular the intensity around 9 = 90 and 270 deg does 
not exhibit a clear trend in the power law creep regime. 
Two experimental observations are nevertheless robust. 
The first one is the systematic concomitance of the con- 
trast C reaching its plateau value with the onset of flow. 
This findings demonstrate that the onset of flow is char- 
acterized by a collective rearrangement of all crystallites 
such that they all orient along the flow direction, hence 
the flow is homogeneous, rather than by wall slippage 
(as observed in the vicinity of the yield stress for a mi- 
crogel paste [IH), or shear localization (as observed in 
adhesive emulsions Q , or wet granular media [l7| ) ■ This 
mechanism is moreover in agreement with the rheology 
experiments (fig. ^ which indicates a bulk mechanism. 
Note in addition that data recorded in tangential ge- 
ometry do not show any significant variation with time: 
isotropic rings are always measured, showing that, in the 
direction perpendicular to the flow, the sample texture 
remains powder-like, and the width of the Bragg peaks 
do not vary, indicating a fixed average crystallites size. 
These observations are consistent with the physical pic- 
ture that only the orientation of crystallites changes in 
the creep and at the onset of flow. This also compares 
favorably with simulations of polycrystals that show ro- 
tation of individual grains under a constant external load 

The second robust experimental observation concerns 
lo, the scattered intensity for 9 = (180±2) deg. I arises 
from the scattering of tubes that are perpendicular to 
the flow (tubes parallel to both the vorticity direction 
and the velocity gradient direction contribute to Iq), and 
resist flow. As shown in fig. [3k, io decreases continuously 
with time and drops down to a smaller and constant value 
once the sample flows. Similar features are systematically 
measured, as shown in fig. [4j where the time evolution 
of Iq for several experiments is reported. In fig. [4] time 
is normalized by i c , and the absolute scattered inten- 
sity by /, the intensity integrated over all angles in the 
flow regime [ijj]. We note that not only the time evolu- 
tion but also the absolute value of Iq are very similar for 
all experiments, which suggests that the evolution of Iq 
reflects a very general mechanism. A peak of Iq is nev- 
ertheless sometimes observed, just before the collective 
reorientation of all crystallites (at t — t c ) and well after 
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a sufficient temporal and spatial resolution would be re- 
quired in order to determine the distribution of waiting 
times between reorientation events and compare it with 
the power law distributions measured in creep fracture 



FIG. 4: Scattered intensity in the direction perpendicular to 
the flow (9 — (180±2) deg) for four distinct experiments, with 
applied stress of 2 (circles and up triangles), 2.5 (squares) and 
3 Pa (diamonds) . The time is normalized by the time t c at 
which the sample flows, and the intensity by the scattered 
intensity integrated over all angles in the flow regime. The 
arrows indicate the normalized fluidization time, Tf/t c - 



the fluidization time. It is presumably due to a transient 
orientation of some crystallites (with tubes perpendicular 
to the flow) required to reach a more favorable orienta- 
tion as the stationary flow begins. 

As pointed out above, the whole time evolution of 
the shear rate (fig. [J) presents strong analogies with 
the strain rate time evolution measured in the exten- 
sion creep of composite heterogeneous materials except 
that, in our case, the acceleration of the rate (fluidization 
regime) is followed by a stationary flow with low viscosity. 
By contrast, in extension creep of heterogeneous materi- 
als, it is followed by the fracture of the sample 1(| 12 1. 
Based on this analogy, we propose a phenomenological 
picture for our experimental observations, which is in- 
spired by the fiber bundle models [2(J, and which ex- 
plains the time evolution of both Iq and the shear rate, 
7. The initial state of our material is powder-like, with a 
random and isotropic orientation of the crystallites. We 
assume that "perpendicular crystallites", i.e. crystallites 
for which the tubes are perpendicular to the flow direc- 
tion (those which contribute to Jo) form some sort of 
percolated network, which resists flow. Upon applica- 
tion of a stress, these crystallites progressively reorient, 
hence Iq decreases, starting from the easiest ones (the 
reorientation being more or less easy, based presumably 
on the local configuration with the neighboring crystal- 
lites). Thus, we expect 7 to decrease as the rate of re- 
orientation decreases, until the network formed by these 
crystallites is not percolated anymore. This would sign 
the onset of the fluidization regime, where 7 increases 
and facilitates the reorientation of the remaining perpen- 
dicular crystallites, hence Iq continues to decrease. Once 
all perpendicular crystallites are reoriented, Iq reaches 
a constant minimum value and the sample is expected 
to flow with a constant 7. To reinforce the analogy be- 
tween creep fracture and onset of flow, experiments with 
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